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1 Introduction 

There exists an impressive list of numerical methods for conformal mapping of polygons W. A 
simple observation allows one to extend this list. Once the angles of the polygon are rational 
multiples of tt, the Christoffel-Schwartz (CS) integral which maps the upper half plane H to the 
polygon is an abelian integral on a compact Riemann surface. The full power of the function theory 
on Riemann surfaces may be applied now to attack the evaluation of the CS integral as well as its 
auxiliary parameters. 

In this paper we consider the simplest case beyond the elliptic one (4 right angles in the polygon) 
described in ffi^ and elaborate this approach to the case of a simply connected polygon with six 
right and one zero angle. Two of the right angles have to be exterior ones (i.e. equal to 37r/2) and 
the vertex with zero angle is infinitely distant. We give the representation for the both mappings, 
heptagon to the half plane and back as certain explicit expressions in terms of genus two Riemann 
theta functions. For the latter there exist a robust and effective method of computation [5] with 
controllable accuracy. Therefore we can guarantee the machine accuracy for the conformal mapping 
uniformly in the polygon/halfplane. 

As usual, there are several auxiliary parameters of the mapping which are determined by the 
geometrical dimensions of the polygon. A good portion of those determining equations are linear 
in this approach. Say, if we map the L-shaped rectangular hexagon, one has to solve just one 
nonlinear equation to determine the mapping. 

This method may be used as a reference for testing the numerical conformal mappings. It 
works equally well for the non-convex polygons, in the presence of spikes of the boundary, narrow 
isthma, boundary elements of different scales, etc. 
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2 Spaces of heptagons 

We consider rectangular heptagon with half- infinite width tt " channel" oriented to the east as it is 
shown in the Fig. [TJ Its sides are either vertical or horizontal. The vertex of the polygon at infinity 
is denoted wq , others are enumerated in increasing order with respect to the natural orientation of 
the boundary (counterclockwise). Two vertexes with the angles 37r/2 are given special names Wa 
and w^, 1 < a < /3 < 6. 

The space Va^ of the heptagons with the given indexes a, l3 is parametrized by the lengths of 
the sides, to which we ascribe signs for technical reasons: 

i^Hs := Ws - Ws+i s = l,2,...,5; (1) 

One obvious restriction on the real values Hi , . . . , He, is the following 

Hi-H3 + H5^TT i^: Im{wi-we)). (2) 

The sign of Hg is negative iff a < s < /3: 

{s + ^-a)is + ^-/3) Hs>0. (3) 

The boundary of the heptagon has no self-intersections which means that the dimensions obey the 
following additional inequalities: 



(a,/3) 


Restriction 


(1,2) 


-H2- 


f iJ4 > 




(1,5) 


-H2- 


f iJ4 > 


when Hi — H^ < 


(2,3) 




f iJs > 




(2,6) 


-H2- 


f iJ4 < 


when - Hs + H5 < 


(4,5) 


Hi~ 


> 




(5,6) 


-H2- 


f iJ4 < 





(4) 



Lemma 1 The heptagons with fixed indexes a and j3 make up a connected space Vap of real 
dimension 4 with the global coordinates Hi, ... ^ He, subjected to restriction the sign rule 
and inequalities Q). g 

Of course, the point from the space Va^p defines a heptagon up to translations in the plane only. 
Those translations may be eliminated when necessary by the normalization e.g. wi :— in. The 
reflection in the real axis induces the mapping Va.p 'Pr-p.r-a which in the above coordinate 
system appears as {Hi,H2, . . . , H^) {Hr^,H4, . . . , Hi). 
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Figure 1: Heptagons from spaces "Pse (left) and P15 (right) 

3 Hyperelliptic curves with six real branch points 

The conformal mapping of the upper half plane to any heptagon from the space Vap may be 
represented by the Christoffel-Schwarz integral. This integral lives on a hyperelliptic curve with 
six real branchpoints. In this section we briefly remind several facts about such curves. 

3.1 Algebraic model 

The double cover of the sphere with six real branch points xi < X2 < ■ ■ ■ < X5 < xe a compact 
genus two Riemann surface X with the equation (of its affine part): 

6 

y^^Y[{x-Xs), {x,y)eC\ (5) 

This curve admits a conformal involution J{x, y) — [x, —y) with six stationary points Ps '.— (xs, 0) 
and an anticonformal involution (reflection) J{x,y) = {x,y). The stationary points set of the 
latter has three components known as real ovals of the curve. Each real oval is an embedded circle 
[2] and doubly covers exactly one of the segments [a;2,a:^3], [a;4,a;5], [xgjXi] 9 00 of the extended 
real line M := K U 00. We denote those ovals as first, second and third respectively. The lift of 
the complimentary set of intervals [xi,X2], [x3,a:4], [0:5, oig] to the surface ([5]) gives us three coreal 
ovals which make up the set of points fixed by another anticonformal involution J J = J J. 

3.2 Homologies, differentials, periods 

We fix a special basis in the 1-homology space of the curve X intrinsically related to the latter. 
The first and second real ovals give us two 1-cycles, ai and 02 respectively. Both cycles are oriented 
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Figure 2: Canonical basis in homologies of the curve X 

(up to simultaneous change of sign) as the boundary of a pair of pants obtained by removing real 
ovals from the surface. Two remaining cycles 61 and 62 are coreal ovals of the curve oriented so 
that the intersection matrix takes the canonical form - see Fig. [21 
The reflection of the surface acts on the introduced basis as follows 

Jas^as, Jbs = -bs, s = 1,2. (6) 

Holomorphic differentials on the curve X take the form 

d-u* = (Ci*x + C2*)y~^dx, (7) 

with constant values Ci*, C2*. The basis of differentials dual to the basis of cycles 

duj := Ssj-, s,j = 1,2, (8) 

determines Riemann period matrix H with the elements 

Usj := [ duj- s,j = 1,2. (9) 



It is a classical result that H is symmetric and has positive definite imaginary part [3] . 
From the symmetry properties ^ of the chosen basic cycles it readily follows that: 

• Normalized differentials are real ones, i.e. Jdus — dus, in other words the coefficients C* in 
the representation ([T]) are real. 

• Period matrix is purely imaginary, therefore we can introduce the symmetric and positive 
definite real matrix i7 := /m(n), 

• Zeroes of the differential du2 (resp. dui) lie on the first (resp. second) real oval. 
3.3 Jacobian and Abel-Jacobi map 

Definition 1 Given a Riemann period matrix H, we define the full rank lattice 

L{\l)^m? + I?= I du, du:= {dui,du2f, (10) 

JHi(X,Z) 

in C'^ and the 4-torus Jac{X) := C^/i(n) known as a Jacobian of the curve X . 



This definition depends on the choice of the basis in the lattice Hi{X,'L), other choices bring us 
to isomorphic tori. 

It is convenient to represent the points u G as a theta characteristic [e, e'], i.e. a couple of 
real 2-vectors (columns) e, e': 

u^\{ne + e'). (11) 

The points of Jacobian Jac{X) in this notation correspond to two vectors with real entries modulo 
2. Second order points of Jacobian are represented as 2 x 2 matrices with Z2 entries. Please note 
that we use a nonstandard notation of theta characteristic as two column vectors written one after 
another. Usually the transposed matrix is used. 

Definition 2 Abel-Jacobi (briefly: AJ) map is a correctly defined mapping from the surface X to 
its Jacobian. 

rP 



u{p) 



du mod L (TV), pi:={xi,0); du := {dui,du2) 



(12) 



From Riemann-Roch formula it easily follows |4] that Abel-Jacobi map is a holomorphic em- 
bedding of the curve into its Jacobian. In Sect. Owe give an explicit equation for the image of the 
genus two curve in its Jacobian. Let us meanwhile compute the images of the branching points 
Ps — {xsj 0), s = 1, . . . , 6 of the curve X: 
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where 11^ and E'^ are the s-th columns of the period and identity matrix respectively. One can 
notice that vector e{u{p)) is constant along the real ovals and e'{u{p)) is constant along the coreal 
ovals. 



3.4 Tiling the Jacobian 

Let us consider sixteen disjoint blocks (tiles) in the Jacobian of a curve filled by the points with 
theta characteristics in the sets 



±1 ±1 
±1 ±1 



/:=(0,1), 
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distinguished by the possible choices of four signs. Jacobian itself is the closure of the union of 
those tiles. 

Chosen an orientation on a real oval of the curve, one can distinguish a component in the set 
x~^M. C X which lies to the left of this oval. We shall denote this disc as IHI+. The surface X with 
all real and coreal ovals removed is a disjoint union of four open 2-discs H+, JH+, JH+, JJH+. 
It turns out that AJ map sends each of those discs to a certain block of the Jacobian and each 
(co)real oval - to a certain 2-torus. This helps us to discriminate points p (z X with the same value 
of the projection x{p) 

Theorem 1 Let the disc IHI+ be chosen in accordance with the orientation of (any of) the a-cycles. 
Then the above four discs are mapped to the following four blocks of the Jacobian: 
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Proof. Symmetries of the normalized abelian differentials with respect to involutions J, J 
guarantee the following equalities: u{Jp) = —u{p); u(Jp) = u{p) since the base point pi of AJ 
map is fixed by both involutions. Therefore we may thoroughly investigate the map on the disc 
H+ only. 

Both components Us{p) of Abel-Jacobi map are CS integrals and send the disc H"*" to the 
rectangles with slots shown in the Fig. [S] Clearly, the real part of Us{p) lies in the interval (0, 1/2) 
which may be reformulated as e'{u{p)) e (/, /)*. 

To study the range of 2-vector e := il^^Im 2u{p) we introduce new differentials dv — {dvi , dv2Y '.= 
—ifl^^du. Those are purely imaginary with normalization J^^ dvs = Sjs- Differential dvi has zeroes 
on coreal oval covering the interval {x5,xg); dv2 has zeroes which project to the interval {xi,X2). 
Again, each abelian integral Vs{p) maps the upper half plane to the rectangle with slot similar to 
that in the left Fig. [3l but now the vertex Vs{xi) lies on the right side of the rectangle, therefore 
e{u{p)) := ~2Re v{p) e (/, /)*. g 

4 Two moduli spaces 

We are going to describe the space of CS integrals corresponding to the rectangular polygons from 
the spaces Vaii- This space is an extension of the underlying space of genus 2 curves with real 
branch points. 
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Figure 3: The images of the upper half plane ui(H"'") (left) and U2(H+) (right) 

4.1 Genus 2 curves with three real ovals 

Each genus two Riemann surface X is automatically hyperelliptic, i.e. it admits a conformal in- 
volution J that fixes six points. This involution is unique. Riemann surface is said to have a 
reflection iff it admits an anticonformal involution J (same surface X may have several anticon- 
formal involutions). Each component of the set of J- fixed points is an embedded circle 2, known 
as a real oval of the reflection. 

Definition 3 Moduli space M2^3 is the space of genus two Riemann surfaces with reflection and 
three enumerated real ovals. Two surfaces are equivalent iff there is a conformal 1 — 1 mapping 
between them commuting with the reflections and respecting the marking of real ovals. 

Let us consider the constructive model of an element X G Ai2^3. 

Necessarily, the involutions J and J of X commute (otherwise J is not unique) and therefore 
J acts on the Riemann sphere X/ J. Once the refiection of the surface acts with fixed points, so 
does the induced reflection of the sphere. For a suitable choice of the global coordinate x on the 
sphere, the reflection of the latter works as a complex conjugation. It is convenient to think of x 
as of the degree two meromorphic function on the surface. Obviously, it maps each real oval to the 
subset of the real equator R := RU 00 of the Riemann sphere. There are three possibilities for the 
image of a real oval: (i) x maps the oval 1 — 1 to the equator; (ii) the oval is mapped 2 — 1 to the 
equator or (iii) an oval is mapped 2 — 1 to the finite segment of it. In the latter case the endpoints 
of the segment are the critical values of the projection x. 

A simple combinatorial arguments show that in the case of three real ovals on X, they are 
projected to three non-overlapping intervals of the real equator of the sphere. The endpoints of 
those intervals are the projections of the stationary points of the involution J of the surface. We 
can give them the unique names in the following way. 
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Projections of the first, second and third oval induce either natural cyclic order of the real 
equator, or the inverse. In the latter case we change the sign of the coordinate x. Now we can 
enumerate the endpoints of the intervals in the natural cyclic order as X2, . . . , a;6 so that the first 
real oval of the surface is a;^^ ([2:2, X3]), the second is x^^{[x4,X5]) and the third is x^^{[xq,xi]). 
Note that an interval of the extended real axis may contain infinity as its interior or boundary 
point. We have shown that our first definition of the moduli space is equivalent to the following 

Definition 4 The moduli space ^^2^3 is the factor of the cyclically ordered sextuples of points 
(a;i,...,X6) from modulo the action of PSL2{M.) (— real projective transformations conserving 
the orientation of the equator) . 

Normalizing the branch points e.g. as 2:4 = 00, X5 = —1, xq = 0, one gets the global coordinate 
system on the moduli space: 

< xi < 2:2 < 2:3 < 00. (13) 

Other normalizations bring us to other global coordinate systems in the same space. Yet another 
global coordinate system in this space is related to the periods of holomorphic differentials. 

Theorem 2 The period mapping ^1{X) is real analytic dijjeomorphism from the moduli space 
AI2K3 to the trihedral cone 

< r2i2 < mm(fiii, r222) (14) 

Proof (sketch). First of all, the matrix ft := Im 11 is a well defined function on the moduli space: 
both choices of the intrinsic basis in the space of integer 1-cycles introduced in Sect. 13.21 bring us 
to the same period matrix. 

The CS integrals Uj{x) := f^^ duj map the upper half-plane 1 — 1 to the rectangles with slots 
shown in the Fig. |31 The dimensions of the rectangles are related to the elements of the period 
matrix where from the inequalities (1141) follow. 

Now we see that all the curves X G A^2]R3 with fixed first column of the periods matrix are 
parametrized by the length I of the slot in the left rectangle of the Fig. [31 In particular, 0,^2 is a 
monotonic function of I. Now one can study the asymptotical behavior of this function and show 
that ^22{l) -> 00 when I [9 and Vt22{l) ^21 when I 1 [TD]. g 

Remark 1 The inverse mapping, from the period matrices to the (suitably normalized) branch- 
points of the curve is also effective. For genus two it is implemented by the Rosenhain formulae 
(see section\^ in terms of theta constants. 
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Remark 2 The cone is strictly smaller than the space of all real positive definite symmetric 
2x2 matrices il. It is known (see e.g. Jl^) that any indecomposable matrix from Siegel genus 
two space is a period matrix of some Riemann surface. The decomposable matrices make up a 
codimension one Humbert variety which is determined by vanishing of at least one of ten even 
theta constants. In our case of real curves the Humbert variety is a real codimension one subvariety 
in the real 3-space and its complement is disconnected. One can check that 

{i^l) = on the edge {^12 — 0} of the cone, 

= on the edge {^12 = ^11} and 

= on the edge {Q,i2 = 5^22}- 
The definitions of the theta constants wiU be given in the Sect. [S] 

4.2 Curves with a marked point on real oval 

The problems of conformal mapping use a slightly more sophisticated moduli space, that of the 
genus two real curves with a marked point on a real oval (see e.g. ^). 

Definition 5 The space of genus two Riemann surfaces with three enumerated real ovals and a 
marked point po ^ Jpo on the third oval we call A^2,iR3- Two surfaces are equivalent iff there is 
a conformal mapping between them commuting with the reflections and respecting the enumeration 
of the ovals as well as the marked point. 

An argument similar to that in the previous subsection shows that we can introduce an equiv- 
alent but more constructive 

Definition 6 By A^2.iII^3 we mean the sets of seven cyclically ordered points (xg, Xi, . . . , Xg) in 
the real equator M. of the Riemann sphere modulo the action of PSL2(M.). 

Here xq means the projection of the marked point to the sphere, other coordinates are the pro- 
jections of the branchpoints of the curve. The natural projection of 7W2.1K3 to the space 7W2K3 
consists in forgetting of the marked point xq. 

Recall that for the element X of the space A^2R3 there is no natural distinction between two 
components of the set x^^H C X until we orient a real oval. The difference arises once we mark a 
point pq 7^ JpQ on a real oval: there is a unique disc H+ C x^^H with pq on its boundary. 

One can introduce several coordinate systems on the space A42.i^3- Fixing three points of 
seven, say a;4 := 00, := —1, xq :— 0, the positions of the remaining four points of the 7-tuple 
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will give us the global coordinate system on the moduli space: 

< xq < xi < a;2 < < cx). (15) 

Other normalizations of the 7-tuple of points result in different coordinate systems on the moduli 
space. It is a good exercise to show that the arising coordinate change is a real analytic 1-1 mapping 
from the 4D cell (fT5| to the appropriate cell. 

Yet another coordinate system on J^2.i^3 is the modification of that related to the periods 
matrix. Three variables flu, ^12, ^22 are inherited from the space ^^2^3 and the fourth is either 

or 1*2, the component of the image of the marked point under AJ map. The integration path 
for u{pq) is the interval of the third real oval avoiding the branch point pg- 

Lemma 2 Each mapping (xq, xi, 2:2, X3) (r^n, rii2, ri22, u^), s = 1,2 is a real analytic diffeo- 
morphism of the cone U5\) to the product of the cone and the interval (0, 1/2) 

Proof. Let us consider the points Xg, s = 0, 1, 2, 3 as the coordinates in the space M2,i^3, other 
branchpoints being fixed. By definition, the period matrix fl is independent of the position of the 
marked point po and the mapping {xi, X2, x^) CI is real analytic diffeomorphism from the cell 
([Tn|) to the cone (HH). Both basic differentials dui, du2 are real and have no zeroes on the third 
real oval containing the marked point. Therefore, both values Us(po) monotonically increase from 
zero to i = i /ai+a2 ^^'^ when the marked point pq moves from pi to pq along the third real oval. 
■ 

The periods of more sophisticated differentials taken instead of dug bring us to yet another 
coordinate systems in the same moduli space of curves with marked points on them. 

4.2.1 ChristofTel-Schwarz differentials 

Let l<a</3<6bea couple of integers labeling the spaces of heptagons. To each element of the 
moduli space A42.i^3 we ascribe the unique abelian differential dwap of the third kind with simple 
poles at the points pq, Jpq, residues resp. —1, +1 at those points and zeros at the branchpoints Pa 
and Pi3 (one of the ways to normalize meromorphic differentials is to ascribe its zeroes at the points 
of nonspecial degree g divisor). The differential will automatically have double zeroes because is 
is odd with respect to the involution J. In the algebraic model ([5]) it takes the form: 

dx 

dWap ^ [x ~ Xa){x - Xjs)— (16) 

where = 11^=1 ~ ^j) ^^-^^ t'^^ 7-tuple oo =: xq, xi, . . . , represents an element of the space 
A^2.iII^3- Christoffel-Schwarz differential ([T5)) can be decomposed into a sum of three elementary 
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ones: 



+ Cidui + C2du2 



(17) 



where dvjp^p^ is a-normalized third kind abeUan differential with poles at Jpo and po; dui and 
du2 are above normalized holomorphic differentials and the constants Ci and C2 are uniquely 
determined by the condition that dwap has zeroes at the branching points Pa and pjj. 

Each element of the heptagon space Vap is the image of the upper half plane under the 
Christoffel-Schwarz map 



Actually, even more is true. 

Theorem 3 Christoffel-Schwarz mapping Wa(j {x) induces a real analytic diffeomorphic map from 
the moduli space A^2,iR3 to the heptagon space Pap- 

Proof. First of all we check that each CS differential dwai3 is real. This implies that the increment 
of the CS map on the boundary of the disc embedded to the Riemann surface is real on the 
real ovals and it is pure imaginary on the coreal ovals. This increment is monotonic between the 
branchpoints but at pq. The image of 9111+ under the CS map is a polygonal line with the same 
sequence of corners as in any heptagon from the space Vap, namely all the corners are equal to 
tt/2 but those at Pa and pp where it is equal to 37r/2. Moreover, this polygonal line has no self- 
intersections. For instance, in the case {a, /3) ~ (2,6) the self-intersection shown in the Fig. H]is 
impossible because the points of the rectangle marked by * have —1 preimages in H according to 
the argument principle. 




(18) 



Figure 4: Self-intersection of the image of real equator under W26 map. 
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We have estabhshed a 1 — 1 correspondence between the moduh space and the space of hep- 
tagons induced by the CS map (|18p . Now wc prove the real analyticity of the direct and the 
inverse mapping. One can easily see that the dimensions of the heptagon are the periods of the 
corresponding CS integral (since CS differential is odd with respect to J): 

^^^3 Ip^ dwafi = iTT + i /^^_^^ dwafi; (19) 

i'^H4, Jl^^^ dWafi = -5 X,^ dwafi; 

*^-^5 Jp^ dwc^p = i J_^^_ dw^p. 
The basic 1-cycles may be separated from the branchpoints of the surface, so the real analyticity 
of the direct map "Pap is clear. It remains to show that the map has the full rank. 

Let us consider the coordinate system in the moduli space such that xq = 00, Xa — xp — 1. 
Assume that the CS-induced map degenerates at a point {xq, xi, . . . ,xq) of the moduli space, then 
there exist a nontrivial tangent vector ^ :— X]i=j^{a ^} Oaf" annihilating all the periods of the 
CS differential dwaf3 at this point. This means the existence of the mcromorphic differential 

l=j^{a,/3} ^ 

with zero cyclic and polar periods on the associated surface and with poles at the branchpoints of 
the surface, except Pa and pp. Then the meromorphic function v{x^ y) := j^^'^^ dv is single valued 
on the surface and it has at most four simple poles along with triple zeroes at the points Pa and 
pp (since dv = 0). This function is identical zero and therefore the above tangent vector is also 
zero. Hence the CS-induces mapping M.2,i^z T^a/s has the full rank everywhere. g 

4.2.2 Auxiliary parameters of the CS map 

Given a heptagon from the space Vap, the corresponding point in the moduli space A^2,iIIi3 may 
be found as (unique as it follows from the Theorem[3|) solution of the following system of equations 
obtained from ([T^ by substituting the decomposition (|17p of Christoffel-Schwarz differential into 
elementary differentials. (We have also used Riemann bilinear relations (|30l) to transform the 
6-periods of elementary differential of the 3-rd kind to Abel-Jacobi image of its poles.) 

2Hi = {Cinn + + 27r(l - 2^°)); 

(20) 

2Hi — — C2; 

2i/5 - -(Cirii2 ^21^22 - ^^A). 
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Here Ci and C2 are real analytic functions on the moduli space defined above as a solution of 
linear system with the nonsingular 2x2 matrix \\du{pa),du{pi3)\\. We see that the auxiliary 
system of equations is " almost linear" with respect to the coordinate system related to the periods 
of holomorphic differentials. It remains to get the effective evaluation of the functions Ci,C2 and 
1*2 as well as the CS integral itself. This is the subject of the next section. 

5 Theta functions on genus two surfaces 

Here we give a short introduction to the theory of Riemann theta functions adapted to genus 
two surfaces. Three problems related to conformal mappings will be effectively solved in terms of 
Riemann theta functions: 

• Localization of the curve inside its Jacobian; 

• Representation of the 2-sheeted projection of the curve to the sphere; 

• Evaluation of the normalized abelian integral of the third kind (which is the essential part 
of the CS-integral). 

Definition 7 Let u G and H G C^^^ be a Riemann matrix, i.e. H = H* and Im H > 0. The 

theta function of those two arguments is the following Fourier series 

6{u,I\) -.^ ^ exp(27riTO*M + Trim'Hm), (21) 

The convergency of this series grounds on the positive determinacy of Im H. The series has 
high convergency rate with well controlled accuracy [5 . Theta function has the following easily 
checked quasi-periodicity properties with respect to the lattice Liji) := lil? + 1?: 

9{u + Hm + m', H) — exp(— iTrm Hm — 

One can say that theta function is a section of a certain line bundle over the Jacobian. In 
particular, its zeroes make up a well defined set known as a theta divisor (9). The Abel-Jacobi 
map transfers the theta function to the Riemann surface where it becomes the multivalued function 
(or a section of a certain line bundle) 

e,{p) :=0(u(p)-e;n), 

which is multiplied by some nonvanishing factors once its argument p goes around the handles of 
the surface. Therefore, its zeroes are well defined. 
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The zero set of the theta function is described by so cahed Riemann vanishing theorems the 
most important of them is the foUowing 4 . 

Theorem 4 (Riemann) Function 0e{p) on a genus two surface either 

(i) vanishes identically on the surface iff e is a certain effectively calculated point K of Jacobian 
(a.k.a. vector of Riemann constants), or 

(ii) has exactly two zeroes qi,q2 such that qi ^ Jq2 and 

e = u{qi) + u{q2) + K mod L{n) (23) 
5.1 Vector of Riemann constants 

We have to determine the vector of Riemann constants K for our particular choice of the initial 
point in the AJ map and the choice of the basis in the homologies. Consider the theta function 
with characteristics which is the slight modification of the above theta. 

Definition 8 

e[2e, 2e']{u, U) := ^ exp(27ri(m + e)*(u + e') + 7:i{m + e)*n(m + e)) (24) 

= exp(i7re*ne + 2i7re*(u + e'))e{u + Ue + e', H), e, e' G M^. 

The matrix argument 11 of theta function is usually omitted when it is clear which matrix we mean. 
Omitted vector argument u is supposed to he zero and the appropriate function of 11 is called the 
theta constant: 

e[e,e'] :=0[e,e'](O,n). 

Remark 3 (i) Theta function with integer characteristics [2e, 2e'] is either even or odd depending 
on the parity of the inner product 4e* • e'. In particular, all odd theta constants are zeroes, 
(ii) Adding matrix with even entries to integer theta characteristics can at most spoil the sign of the 
theta function. Hence, the binary arithmetic plays a great role in the calculus of theta functions. 

Lemma 3 

K = u{pz) + u{pz) (25) 

Proof. We see from the table is subsection 13.31 that the points and correspond to the 
odd characteristics. Remark [3I^i) implies that the function 6{u{p)) will have zeroes at the points 
p = P3,P5 on the surface. This function cannot be identical zero on the surface: otherwise from 
part (i) of Riemann vanishing theorem it would follow that K ~ Q. Now one can check that the 

14 



function 9e{p) with the shift e := u{p^) + u(p^) + K is not identical zero (since ^ Jpb) and 
vanishes at three points p = pi,p3,p5. 

An alternative (ii) of Riemann vanishing theorem suggests that K ~ ^(ps) + ^(ps) mod L(Il). 
In other words, K corresponds to the odd characteristics 

Remark 4 It is convenient to represent integer theta characteristics as the sums of AJ images of 
the branchpoints, keeping only the indexes of those points: [sk..l] means the sum modulo 2 of the 
theta characteristics of points u{ps), u{pk), ■ ■ ■ ,u{pi), e.g. [35] is the vector of Riemann constants 
represented by theta characteristics. 

5.2 Image of Abel-Jacobi map 

The location of genus 2 curve embedded to its Jacobian may be reconstructed by solving a single 
equation. 

Theorem 5 (Riemann) A point e of Jacobian lies in the image u{X) of Abel-Jacobi map if and 
only if 

e\ '° j(e) = 



Proof. 1. The function 0{u{p) — K) vanishes identically on the curve and so does the function 

om{<p))- 

2. Conversely, suppose that 0[35](e) = or equivalently, 9{e + K) — 0. The function 6e'{p) with 
the shift e' :— e + K either vanishes at two points pi, p' and in this case e = u{p'), or identically, 
then e = = u{pi). g 

5.3 Projection to the sphere 

Any meromorphic function on the curve my be effectively calculated via the Riemann theta func- 
tions once we know its divisor. Take for instance the degree 2 function x on the hyperelliptic curve 
([5]). This projection is unique if normalized e.g. as follows: x(ps) — 0, x{pj) = 1, x(j)i) = oo, s,j,l 
is a positive triple from the index set {1,2,..., 6}. 

With the use of the transformation rules one immediately checks that the following function 
is single valued on the curve: 

_ 9^[sk35]{u{p),n) 
"^^P' - 0mk35]{u{p),n)' ^ ' 
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Now with the help of Riemann theorem one checks that the numerator of the function has double 
zeros at the points ps, Pk while the denominator has double zeroes at pi, pk- Therefore, x{p) 
differs by a constant factor from the above normalized projection: 

fl2[/fcj35] 9^[sk35]{u{p)) 



x{p) = ±- 



k ^ s,l,j, 



(26) 



"6i2[sfcj35] 0^[lk3b]{u{p)) ' 
where the sign ± in the latter formula depends on the parity (+ even /— odd) of the scalar product 
^U) • (^'(s) + ^'(0)j where ne(s) + e'(s) :— u{ps). In the case k — j the normalization have to use 
the L'Hospital rule and we get a slightly different factor. 

Also, the standard normalization x{pa) — —^]x{pp) — l;x{po) = oo brings us to a little more 
awkward expression: 

«2[/3j35](wO) e^[aj35]{u{p)) 



x{p) 



±2- 



1, 



3 + OL, 13, 



02[a/3j35] n±^b'35](u(p)±«o) 
again, the sign in front of the fraction depends on the parity of the scalar product e(/3) ■ e'(a), 
where e(s), e'(s) is the representation of the half-period u{ps) by integer theta characteristic. 



5.4 Rosenhain formulae; Humbert variety 

Putting in the above formula (j26l) the half-periods corresponding to the branchpoints, we get the 
effectively computed expressions for the latter in terms of the period matrix. For instance, if we 
normalize the projection x{p) as x{pi) — 0, x{p2) — 1, x{p^^ — oo, we get the following expressions 
for the remaining branch points in terms of theta constants: 



2^3 



02 


00 


02 


00 




00 




01 


02 


01 


02 


01 




00 




01 



2^5 



02 


00 


02 


10 




00 




10 


02 


11 


02 


01 




11 




01 



(27) 



Choosing other normalization for the projection, we get certain expressions for the cross-ratios of 
the same set of branch points, e.g. 



\- Xz 



02 


10 


02 


10 




00 




01 


02 


01 


02 


01 




00 




01 



1 — X4 = — - 



02 


00 


02 


10 




10 




01 


02 


01 


02 


11 




00 




11 



1 - = - - 



02 


00 


02 


10 






10 




00 




02 


11 


02 


01 






11 




01 





(28) 



Comparing the formulae for the same branch points from the equations ((27|) and ((28)) we get certain 
relations for the theta constants which are the consequences from the Riemann theta identities. 
The expressions for the branch points of genus two (hence hyperelliptic) curves in terms of theta 
constants appeared in Rosenhain's work |11| . 

Theorem 6 /IT/ For nonsingular genus 2 curve all 10 even theta constants 0[e, e'], e* • e' G 2Z, 
are distinct from zero. 
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Proof. Choose three distmct branchmg pomts p^, pj, pi on the curve. Since ps ^ Jpj, the 
function Oe{p) does not vanish at p ~ pi when e = u{ps) + u{pj) + K. The value Oe{pi) vanishes 
simultaneously with the theta constant 6[sjl35\. Once the indexes s, j, I vary in the set {1, 2, . . . , 6}, 
the thcta characteristic [sj735] runs through all 10 even characteristics, g 

Humbert variety in complex Siegel space is the locus where at least one of even theta constants 
vanishes. This is exactly the locus of Riemann matrices that are not period matrices of any genus 
two Riemann surface. 

5.5 Third kind abelian itegral 

On any Riemann surface there exist a unique abelian differential of the third kind dvrq with simple 
poles at two prescribed points r,q only, residues +1, —1 respectively and trivial periods along all 
a— cycles. It has certain physical meaning in terms of flow of inviscid incompressible fluid on the 
surface with a source at r and sink at q. Our approach to effective conformal mapping grounds on 
the fact that the integral of this differential can be expressed in closed form. 

Theorem 7 (Riemann) Choose any point s G X , then for any two points r,q ^ Js of the surface 
holds the representation: 

r, . e(e + u(p) - u(r)) , e[e,e'](u{p)~u(r)) 

vr,{p) := / dvr, = log Y' ^ = log ' , Y[ H( + 29 

J* d{e + u{p) - u{q)) e[e,e'\[u{p) - u{q)) 

where theta characteristic [e, e'] corresponds to the zero e :— —u{s) — K of theta function (say, it 
may be any odd integer one). 

Proof. Consider the following function of variable p £ X: 

6{e + u{p) -u{r)) 
e{e + u[p) - u[q)) 

One checks that it is locally holomorphic on the surface: the theta function in numerator vanishes 
at the points r, s, the denominator vanishes at q, s. Moreover, this function is single valued on the 
surfacqj. To prove this we use the Riemann bilinear identity 

dVrq = 2TTi / duj (30) 

and the transformation properties (j22p of theta functions. Therefore, it is a constant independent 
of p. ■ 



^provided the branches for u{r), u{q) in II29I I are properly chosen: the integration path from r to the initial point 
pi and further to q may be deformed to the one disjoint from a— and b— cycles 
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Example: CS integral from the section l4.2.1l mav be represented as follows: 

, , , e[3]{u{p) + u°) ^ , , ^ , , 

W^p{p) = ^Og g[3](^(p)_^0) + ClUlip) + C2U2{p), 

where u{p) := {ui{p) , U2{p)Y ', := u{po) and the constants Ci,C2 are obtained from the system 
of linear equations dwapiji-y) — 0, when 7 = a, /3. 



6 Algorithm of conformal mapping 

Based on the formulae of the previous section, we can propose the algorithm of conformal mapping 
of the heptagon to the half-plane and vice versa. First of all, given the heptagon we have to 
determine the corresponding point of the moduli space 7^2,11X3. 

6.1 Auxiliary parameters 

Given (related to side lengths) coordinates Hs of the heptagon, we have to determine seven real 
parameters: the (imaginary part of period) matrix i7, the image :— u{po) of the marked point 
Po in the Jacobian of the curve and real 2- vector C := (Ci, C2). Those parameters give a solution 
to a system of seven real equations: 

d0[35](u,*17)Ad(log^}|^^^±447T + ^-") = O' when u^u{p^),u{pf,) (31) 
0[S\[u ~ u'^ 

which means that CS differential dwa.p has zeroes at the points Pa and pp; 

e[ib]{u°,in) ^Q, (32) 

which means that the point vP lies on the AJ image of the curve in the Jacobian and finally 

Hi = i(CiOii + 02^12 + 27r(l - 2u\)); 

^ (33) 
H4 = — 5C2; 

H5 = -^{Cirtu + €2^22 - 47ru§). 
which specify the side lengths of the heptagon. 

Lemma 4 Let the side lengths Hi , H2 . . . , H^ satisfy the restrictions described in lemma [71 then 
the system of seven equations I131\) . II 3 i33\) has a unique solution fl, vf^ , C £ in the domain 
determined by inequalities: 

< rii2 < min(f7ii, $122), 
0<M°<l/2, s = l,2. 
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Proof. The existence and the uniqueness of the solution to ([5T|) - ([55]) in the specified domain fol- 
lows from the existence and the uniqueness (up to real projective transformations) of the conformal 
mapping of a given heptagon to the upper half plane. g 

Remark 5 Essentially we have just three nonlinear equations to solve which can be done by Newton 
method with parametric continuation. Indeed, take any point (f2, Uj*) in the moduli space, then 
solving just one nonlinear equation iS^) we get the whole vector vP , then solving two linear equations 
\31]) we get the constants C and substituting the values fi, , C to the last four equations 13!^) 
we get the dimensions of the appropriate heptagon. So we got the correspondence of two points: 
one in the space Va^p and another - in the space A^2,iIR3- Now, any given point in the heptagon 
space may be connected to the obtained one by a path of at most two linear segments. Those may 
be effectively lifted to the moduli space by the Newton method. 

6.2 Mapping heptagon to half plane 

Suppose a point w* lies in the normalized heptagon (i.e. with the vertex wi = in). Consider a 
system of two equations 

8[3\{u' — u" ,iflj ' (34) 

9[35]iu*,in) = 0, 

with respect to the unknown complex 2-vector u*; real symmetric matrix fl and real 2- vectors u^, 
C being the solutions of the auxiliary system ([5T|) - 

Lemma 5 System fi34\ l has the unique solution u* with theta characteristic from the block 

Proof. The left hand side of the first equation in (IM]) restricted to the AJ image of the marked 
surface {X,pq) in represents the abelian integral Wap. The abelian integral is not a single valued 
function on the surface X . However there is a single valued branch of the integral in the disc H+ 
which maps the latter 1 — 1 to the given normalized heptagon. Therefore, the point w* has the 
unique preimage in the disc C X. Abel-Jacobi image u* of this point satisfies two equations 
([M)) and lies in the specified block of the space C^. g 

Remark 6 It is clear from the reflection principle for the conformal mappings that the set of two 
equations \34-^ may have many solutions. We use theta characteristic to single out the unique one. 

Substituting the solution u* to the right hand side of the expression ([^5]) we get the evaluation 
at the point w* of the conformal mapping x{w) of the heptagon to the half plane with normalization 
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Ws, Wj, Wl — >■ 0, 1, 00. 



e^[skj35] 0^[lk35]{u*) ' 



k ^ s,l,j. 



(35) 



6.3 Mapping half plane to heptagon 

Conversely, given a point x* in the upper half plane we solve a system of two equations with respect 
to a complex 2-vector u* with characteristics from 



-I I 
-I I 



0^[lkj35] 9^lsk35]iu*) 
e-'lskj35\ 6i^[;fc351(«*) 



6»[35](w*) = 



(36) 



and substitute this solution to the formula 



log rrf + C ■ u := w 



to get the image of the point x* in the heptagon. 



(37) 



7 Conclusion 

It often happens in classical mechanics that the problem becomes much simpler when considered 

in adcqiiatc coordinate system. We sec that the same is true for the conformal mapping: tak- 
ing elements of the period matrix as independent coordinates on the moduli space of ciirvcs, we 
have simplified the system of equations for the auxiliary parameters of conformal mapping. The 
mappings themselves are now closed analytic expressions involving effectively computed Riemann 
theta functions. 

Our approach to conformal mapping of (rectangular) polygons can readily be generalized to 
other types of polygons. Of course, for involved polygons the the genus of the associated curve 
becomes higher which makes the algorithm more complicated. For the set of auxiliary parameters 
one has to solve the (hyperelliptic) Schottky problem additionally. As well, more efforts have to 
be spent to localize the AJ image of the curve inside the multidimensional Jacobian. 



8 Appendix: asymptotic of periods 
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